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Abstract 

A similarity transformation is obtained between general population matri¬ 
ces models of the Usher or Lefkovitch types and a simpler model, the pseudo- 
Leslie model. The pseudo Leslie model is a matrix that can be decomposed 
in a row matrix, which is not necessarily non-negative and a subdiagonal 
positive matrix. This technique has computational advantages, since the so¬ 
lutions of the iterative problem using Leslie matrices are readily obtained . 
In the case of two age structured population models, one Lefkovitch and an¬ 
other Leslie, the Kolmogorov-Sinai entropies are different, despite the same 
growth ratio of both models. We prove that Markov matrices associated to 
similar population matrices are similar. 
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1. Introduction 

This article deals with classic discrete structured models for linear popula¬ 
tion dynamics j3,[^ such as Leslie matrices and Lefkovitch or Usher matrices. 
Giving A, a non negative n x n matrix and a population vector which 
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components are the fractions of the population at each age or stage, the dy¬ 
namical system that gives the population vector at any positive time k + 1 
is given by 

Xfc+i = v4xfc, with initial condition xq. 

Obviously the solution is given by the powers of A 

Xfc = A^Xq. 


In this paper we prove that there is a similarity transform that converts 
the complicated dynamics of the so called Usher or Lefkovitch matrices to the 
simpler study of matrices which are Leslie matrices or pseudo-Leslie matrices, 
a concept that we introduce in this paper. 

The paper is organized in three sections, in the second we introduce 
pseudo-Leslie matrices and prove the main theorem. In the third section 
we present some consequences of interest in population dynamics, namely on 
the similarity of Markov matrices associated to similar population dynam¬ 
ics matrices and obtain transformation rules for corresponding stationary 
distributions. 

2. Main theorem 

In age structured population dynamics one divides the population in 
classes [1, [^. When we consider size classes or stage classes instead of pure 
age classes we have a structured population model with dynamics given by 
the linear equation 

^ 71+1 ( 1 ) 

where x„ is a non negative structured absolute population vector, or a pro¬ 
portion of individuals in each class and £ is a matrix such that 
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usually called Usher (in the classic reference [^) or Lefkovitch matrix in [^. 
The coefficient fj is called the fertility rate of class j > 1, the coefficient 
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bk > 0, for any k = — 1, is the transition rate from class A; — 1 to 

class k and the q the rate of individuals that remain in class 1. Along this 
paper we assume that /„ > 0, assuring that C is irreducible j^. 

The coefficient /i can be decomposed in /i + cq, i.e., a fertility rate and a 
permanency rate. Since this decomposition has no influence on the similarity 
transformation we do not split /i. One must keep in mind the biological 
meaning of this coefficient. 

The solution of the problem is given by the powers of £, given the non¬ 
negative initial condition xq 


x„ = /I’^Xq. 

A Leslie matrix is a matrix of the type 
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where all the entries (pj are non-negative and all bj are strictly positive. The 
Leslie matrix can be decomposed in two matrices 

L = R + B, 

where 
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When the entries <f)n of the hrst row of R are real numbers, not restricted to 
the non-negative case, we say that L is a pseudo-Leslie matrix. Obviously this 
class of matrix does not have an immediate biological correspondence when 
some of its entries are negative. That poses no problem in the framework of 
this article, since L is merely used as a computational instrument. 
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To state the main theorem we dehne the sums of products of p factors T?, 
where i = 1 , ...,n denotes the row index of a given n x n Lefkovitch matrix 
C 
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if p = 0 
if n — z < p 


For the products of the transition rates bi,6n-i of £ we use the notation 


M = 


3 

if z < j < n — 1 

k=i 

1 if j = z — 1 


Now we introduce an upper triangular matrix S and a pseudo-Leslie matrix 
L defined by 
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We are now in position to state the main result of this work. 

Theorem 1. For any Lefkovitch matrix, C, one has S~^CS = L where S 
and L are the matrices defined above. 


The following lemma is used in the proof of theorem [T] 

Lemma 2. If C is a nxn Lefkovitch matrix, then Tf^^ = Q-iTf + for 
all p > 0 and n > i > 1. 

Proof. As TO = 1, T] - = q.i and Tf = = Tf+i = 0 for p > n - z, 

the proof is obvious for p = 0 or p > n — z . So we may assume 0 < p < rz — z. 
If 0 < p < rz — z, then 


rp; = (- 1 )'’+* 

n—l>ip+i>--->i2 >*!>*—1 

= (-ir‘ E 


CjiCjj • • • Cjp+i 


Q 1 Q 2 ■ ■ ■ ^ip+i 


n—l>ip+i>--->fi2>n>j 

- (- 1 )' E 


Ci—lCi2 ■ ■ ■ Cip_|_i 

n—l>ip+i>--->i 2 >i 

= rp' - c_.rf. 

Finally, assume that 0 < p = zz — z. In this case, as Ff"'"^ = 0, one gets 

c.-irf + rp; = 


Cj_i ( 1) CjCj+1 . . . Cn—± T ( 1) Cj—iCj . . . Cji—i 
= 0 = Ff+\ □ 


We are now in position to prove the main result. 

Proof of theorem [H In order to prove the equality jCS = SL, we begin by 
computing SL. As Si^i = 1 and Sij = 0 for z > j, one has 




Si,i(f)n if j = rz 

T j ^ rz 

4>n if z = 1, j = rz 




= b 


Sij+iOj 

0 


if z = 1, j < rz 
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if rz > j > z > 1 
otherwise 
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As = 0 and A{ = A{ ^bj one has 



and 





finally we get 


pj+l-i pi+i-* 

Si^j+ihj = bj = , for n> j >i. 


Thus, we may write 




(Si)., = <! 


k=l 
6 , 

pj+i-i 

_i_ 

0 


if f = j + 1 
an>j>i>l 
otherwise 


Notice that since T”’'"^ * = 0 for all i, we finally arrive at 
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otherwise 



Next we compute CS. As Si^i = 1 and Sij = 0 for i > j, one has 
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1 
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(3) 


Now, using lemma [2] we see that ([2]) and ([2D are the same, which completes 
the proof. □ 

The dynamical system o can be solved using the easily computable 
powers of L 

x„ = £^xo = S-^L^Sxo. 

Since £ and L are similar, they share the same spectrum and the Perron- 
Frobenius Theorem still holds for L in what concerns the existence of a simple 
dominant positive eigenvalue. Using a generating function and formal power 
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series obtained in [l(| or the classic Jordan canonical form, it is always possible 
to obtain the powers of L. The eigenvectors of £, will be studied in the next 
section. 


3. Sinai Kolmogorov entropy, Markov matrices and stationary dis- 
tribntions 

In this sectiom using a simple example, we show that the Kolmogorov- 
Sinai entropy [3|, y, y, 1^ is not an algebraic invariant. We also establish that 
two Markov matrices associated to population dynamics similar matrice^ 
are similar. Finally, we establish a transformation rule for the two station¬ 
ary distributions of Markov matrices associated with two similar population 
matrices. 

Given two matrices, one of Lefkovitch type and the other of Leslie typ^, 
with the same growth rate, they can have different Sinai-Kolmogorov en¬ 
tropies as we see in the following example. 


Example 3. Let 


C 


we have the similarity matrix 


1 3 

0.4 0.55 ’ 


S 


1 -1.375 
0 1 


and a Leslie matrix L similar to L, which is 

_ r 1.55 1.625 ■ 

^ [ 0.4 0 

The Perron-Frobenius dominant eigenvalue is X = 1.89331 both for L and 
C. The Markov matrix corresponding to a population matrix A is 

obtained using the relations 



XUi 


^Under very general conditions. 

^We consider a true non-negative Leslie matrix to establish this conclusion. 










where A is the dominant eigenvalue of A, and the column vector u = ^ > 

0 is the Perron-Frobenius right eigenvector of A. (The left eigenvector will 
he called the line vector v = ^). For the Lefkovitch matrix C we get 

the associated Markov matrix 


pC 


0.528175 0.471825 
0.709504 0.290496 ’ 


the stationary distribution of P^ is = \_ 0.600598 0.399402 ]. The pop¬ 
ulation Sinai-Kolmogorov entropy f^] is 


2 

He = -Y^ logpg, 


where pfj are the entries of P^ and Trf are the components of the station¬ 
ary distribution vr^ of P^ (the left eigenvector associated with the Perron- 
Frohenius eigenvalue 1 of P^, such that tt^P^ = n^). Doing the same com¬ 
putation for L we have 


2 

Hl = 


where P^ is the matrix with entries ph, the Markov matrix associated to L 
is 


L _ 0.818671 0.181329 

^ [ 1 0 

The stationary distribution of P^ is = [ 0.846504 0.153496 ] and the 
entropies of C and L are different, respectively He = 0.656027 and = 


0.400738. 


The Markov matrices P^ and P^ associated to L and C are also similar, 
with the same eigenvalues as we will see below. This result can be stated 
in the general context of similar matrice^ under the following hypothesis, 
which are assumed until the end of the paper: 

1. £ is non-negative and irreducible, therefore has the dominant eigen¬ 
value A, and associated left and right positive eigenvectors t and w, 
respectively. 


®Not necessarily Lefkovitch, Usher or Leslie matrices. 
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2. L and £ are similar, related by the invertible similarity matrix S, snch 
that CS = SL. 

3. L, not necessarily non-negative, has right and left eigenvectors, respec¬ 
tively u and V, associated to A with all entries positive. 

The right eigenvector of L associated to the dominant eigenvalne A 

Lu = Au 


is related to the right eigenvector w of £ by the transformation rnle w =5'u, 
since 

£S'u =AS'u y^£w =Aw. 

The same happens for the left eigenvector v of £ 

vL = Av 


and the left eigenvector t = \S ^ of £, since 

v£-^£=Av£-^y» t£=At. 

The Markov matrix associated with L is given by its entries 

L ^ij '^j 

On the other hand, the Markov matrix associated with £ is given by 

Pi = 


r _ ^ij U)j 


XWi 


The stationary distribntion of is 


TT^ = 


[ ViUi V2U2 . . . VnUn ] 


vu 


where vu is a compact notation for the inner prodnct of the line vector v 
and the colnmn vector u. The stationary distribntion of £ is 


TT^ = 


[ tiWi t2W2 ■ ■ ■ tnWn ] 


tw 


It is possible to prove that the Markov matrices and P^ are similar. 
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Proposition 4. and are similar if L and L are similar. 
Proof. One defines the sqnare matrices U and W snch that 


u = 

Ml 

U2 


, w = 

Wi 

W2 







Wn 


with all Wj 7 ^ 0 and Wi ^ 0, the inverses of P and W are 


- J_ 


■ J_ 

Ul 


W\ 

J_ 



U2 

, = 

W2 



1 

- 


yon _ 


With this notation consider the transformations 

P^ = \u-^LU and P^ = ^W-^CW, 

A A 

where A 7 ^ 0. 

Now, it is straightforward to prove that P^ and P^ are similar 
P^ = ^w-^CW = \w-^SLS-^W. 

A A 

On the other hand 

P^ = \u-^LU. 

A 

Therefore, XQ and XP are similar, since both are similar to L. Explicitly 
L = XS-^WP^^W-^S = XUP^U-^ 


or 

pL = U-^S-^WP^W-^SU, (4) 

as desired. □ 

We can prove that vr^ is a stationary distribution of P^ using matrix 
notation. 

Proposition 5. The row vector is a stationary distribution of P^. 
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Proof. Using the left eigenvector t = [ U t 2 ... ] of £, we define 

diagonal matrix 

T = 


We have 


CpC ^ -r ^ ^ _ _ ^1 TWW-^CW 

Atw L J 

[ ti t2 ... tn ] £fU, 


Atw 


since t is a left eigenvector of £ we have 


^2 ... tn ] W 

TW 


tw 

= TT^. □ 


Using analogous techniques we obtain the relation between the two sta¬ 
tionary distributions of and P^. 

Proposition 6. The stationary distributions vr^ and tt^ are related by 

= n^W-^SU 


Proof. From (jl]) we have 

= z-^p^z, 

where Z = W~^SU. In that case the stationary distribution tt^ is given by 
the relationship 

pL L 

Tl r = TT , 

SO 

-I^Lz-lpCz = T^L^-lpC ^ 

which means that 

C L y — 1 

TT = TT Z , 

as desired. □ 
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Remark 7. All the results in this section apply to the case of an irreducible 
Lefkovitch matrix C and a similar pseudo-Leslie matrix L, since any matrix 
of the form 


01 

02 

4^71—1 

0n 

bi 
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0 

0 

0 

b2 

0 

0 

0 

0 

1 
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with positive coefficients bj and with the dominant eigenvalue X has the pos¬ 
itive right eigenvector 


r A? 1 


r 1 1 



bi 



60)2 

A2 

— 

A2 



bib2---b„-i 

L A"-1 J 


L A"-i J 


The similar Lefkovitch matrix 


fi f2 /s 

bi Cl 0 

0 b2 C2 


n—1 Jn 
0 0 

0 0 


0 0 0 • • ■ Cn-2 0 

0 0 0 • * • bfi—\ Cyi—\ 


is always irreducible if fn > 0 and all the bj are positive, Therefore, 
similar Lefkovitch and pseudo-Leslie matrices, C and L, satisfy conditions\J\, 
d and\^ 
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